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SUMMARY 

A magneto -hydrodynamic model of boundary layers at the Core-Mantle Boundary (CMB) 
is derived and used to compute the viscous and electromagnetic torques generated by the 
Earth's nutation forcing. The predicted electromagnetic torque alone cannot account for 
the dissipation estimated from the observations of the free core nutation. The presence of 
a viscous boundary layer in the electromagnetic skin layer at the CMB, with its additional 
dissipative torques, may explain the geodetic data. An apparent Ekman number at the top 
of the core between 2 and 4 1CT 11 is inferred depending on the electrical conductivity of 
the mantle. 

1 INTRODUCTION 

Detailed models of coupling at the Core Mantle Boundary (CMB) have been put forward to explain 
the more and more accurate measurements of the nutations of the Earth (©warning Citation 'wahr81' 
on page 1 undefined; ©warning Citation 'deha97' on page 1 undefined; ©warning Citation 'math02' 
on page 1 undefined). The nutations of the Earth induce a differential rotation, about an equatorial 
axis, between the mantle and the core (©warning Citation 'sasa80' on page 1 undefined; ©warning 
Citation 'buff92' on page 1 undefined). This differential rotation at the CMB generates both a viscous 
torque (©warning Citation 'gree68' on page 1 undefined; ©warning Citation 'lope75' on page 1 
undefined; ©warning Citation 'roch76' on page 1 undefined) and an electromagnetic torque due to the 
shear of the poloidal magnetic field lines (©warning Citation 'roch60' on page 1 undefined; ©warning 
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Citation 'toom74' on page 1 undefined; ©warning Citation 'sasa77' on page 1 undefined). Buffett and 
his colleagues developed sophisticated models of the electromagnetic torque at the CMB (©warning 
Citation 'buff92' on page 1 undefined; ©warning Citation 'buff93' on page 1 undefined; ©warning 
Citation 'buff02' on page 1 undefined) in order to fit the spatial geodetic observations. First, Buffett 
(©warning Citation 'buff92' on page 1 undefined) introduced a weak magnetic field theory where 
the Lorentz forces associated to the skin magnetic effect are too small to generate any motion in the 
boundary layer. His magnetic analysis requires the presence of a very good electrically conducting 
layer in the lowermost mantle (same electrical conductivity as the core) to get an adequate amplitude 
of the torque. Moreover, Buffett (©warning Citation 'buff92' on page 1 undefined) invoked a enhanced 
magnetic field at the CMB (4 times larger than the observed one) to account for the small scales of 
the magnetic field. The value of the small scales of the magnetic field at the CMB (spherical harmonic 
degree I > 13) cannot be measured at the surface of the Earth because the crustal magnetic field is 
dominant at these wavelengths (©warning Citation 'blox95' on page 1 undefined; ©warning Citation 
'stac92' on page 1 undefined). He estimated their effect using an extrapolation of the low-degree non 
dipole part of the poloidal magnetic spectrum to higher degrees. Then, Buffett (©warning Citation 
'buff93' on page 1 undefined) investigated the role of a toroidal magnetic field on the electromagnetic 
torque at the CMB. Its effects are weak and do not increase the dissipation of magnetic origin at the 
CMB. Moreover, his results are rather speculative as measurements of the toroidal magnetic field in 
the Earth's core are not available. Buffett et al. (©warning Citation 'buff02' on page 1 undefined) 
improved the 1992's model by relaxing the weak field approximation. Thus, they solved the inviscid 
dynamics of the skin layer in the presence of Lorentz forces. The ratio of the velocity induced by the 
Lorentz forces in the skin layer and the velocity jump at the CMB is of the order of the Elsasser number 
(defined below). Its value, at the top of the core, is comprised between 0.1 and 1 so that the weak 
field approximation is not valid. The presence of this dynamical effect reduces the amplitude of the 
electromagnetic torque at the CMB. This is the reason why, in order to fit the improved observational 
constraints (©warning Citation 'math02' on page 1 undefined), Buffett et al. (©warning Citation 
'buff02' on page 1 undefined) invoked the presence of a constant magnetic field modeling the non 
dipole component (small scale magnetic field), three times greater than the dipole value at the CMB. 
In all their studies, Buffett and his colleagues introduced a thin electrically conducting layer at the 
base of the mantle. Its presence remains necessary to get the correct amplitude of the electromagnetic 
torque. 

For rapidly rotating fluids, viscosity plays a role mainly in thin boundary layers, the so-called 
Ekman layers (©warning Citation 'gree68' on page 2 undefined). The depth of these layers is \Jv/VL 
where v is the kinematic viscosity of the fluid and the angular velocity of the Earth. As the magnetic 
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parameter 


value 


R 


core radius 


3.48 10 6 m 


fl 


rotation rate of the Earth 


7.29 1CT 5 rad s^ 1 


o 

r 


density 


10 4 kg m~ 3 


n 


magnetic diffusivity of the core 


1.6 m 2 s _1 


TIM 


magnetic diffusivity of the mantle 


1.6 - 1600 m 2 s" 1 


V 


kinematic viscosity of the core 


7.0 10~ 6 m 2 s^ 1 


Bo 


magnetic field at the CMB 


0.46 10" 3 T 




coupling constant at the CMB 


-1.85 10~ 5 


E 


Ekman number 


8.0 10" 15 


E m 


magnetic Ekman number of the core 


1.8 10" 9 


E M 


magnetic Ekman number of the mantle 


1.8 10~ 9 - 1.8 10" 6 


A 


Elsasser number ~^ 


0.14 


p 

1 m 


magnetic Prandtl number - 


4.5 10~ 6 



Table 1. Physical properties and associated dimensionless numbers used in this study. 



skin depth is where r\ is the magnetic diffusivity of the core, the ratio of the two lengths is 

given by \JP m where P m = v/rj is the magnetic Prandtl number. Table Q] contains the values of the 
molecular diffusivities for the core (©warning Citation 'poir94' on page 2 undefined). We evaluate 
P m = 4 10~ 6 in the core, making the viscous layer 500 times thinner than the magnetic skin layer. 

Recent numerical simulations of the geodynamo have been successful in reproducing some fea- 
tures of the magnetic field of the Earth (©warning Citation 'dormOO' on page 3 undefined). They have 
in common to use a very high viscosity (Ekman number greater than 10 -6 ) so as to avoid numerical 
resolution problems. Glatzmaier & Roberts (©warning Citation 'glat95' on page 3 undefined) advo- 
cated the use of an eddy viscosity for dynamical core modeling as it is generally done in numerical 
modeling of the oceanic or atmospheric sciences (©warning Citation 'pedl87' on page 3 undefined). 
Brito at al. (©warning Citation 'brit04' on page 3 undefined) have found evidence of apparent vis- 
cosity from an experiment of thermal convection in a rapidly rotating spherical shell filled with water 
using a spin-up technique. They interpret their observations by arguing that turbulent motions in the 
bulk of the core increase the efficiency of the exchange of angular momentum between the Ekman 
layers and the geostrophic volume. These non-linear effects at small scale may be modeled by an 
eddy viscosity at large scale. A turbulent viscosity at the top of the core between 10 -4 m 2 s _1 and 
10 _1 m? s^ 1 is possible. Such eddy viscosities increase the magnetic Prandtl number and decreases 
the ratio between the "viscous" and magnetic layer depths. Under these conditions, viscous effects 
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have to be incorporated in the dynamical equation of the layer. With such a theory, the quality factor 
of the free core nutation, deduced from the geodetic data, is a constraint on the apparent viscosity et 
the top of the core for the diurnal frequency. 

Recently, two related studies (©warning Citation 'mafh05' on page 4 undefined; ©warning Cita- 
tion 'palm05' on page 4 undefined) have been published. Their approaches are very similar to the work 
presented here and lead also to the prediction of a viscosity value at the top of the core from nutations 
data. From their own data analysis, Palmer & Smylie (©warning Citation 'palm05' on page 4 unde- 
fined) use an approximate viscous model to infer a viscosity. Mathews & Guo (©warning Citation 
'math05' on page 4 undefined) introduce a magneto- viscous model similar to ours and determine the 
viscosity from the observational data analysis of Mathews et al. (©warning Citation 'mafh02' on page 
4 undefined). Both papers give a value of viscosity which is close to the one proposed in this paper. 
However, our analysis proposes a complete calculation with all spectral components of the magnetic 
field with different extrapolated tendencies for the hidden part (I > 13) of the magnetic field at the 
CMB. Moreover, a physical description of the magnetic and viscous boundary layers is shown. We 
also give a complete study of the variations of the electrical conductivity at the base of the mantle 
which enables us to invert the observational data to obtain trade-offs between the mantle electrical 
conductivity effect and the viscous effect at the CMB. 

This paper presents a derivation of a magneto-hydrodynamic boundary layer attached to the mantle 
taking into account the Lorentz, Coriolis and viscous forces (section 2). In section 3, we discuss 
the influence of the geometry (small scales) and amplitude of the magnetic field at the CMB on the 
electromagnetic torque. Section 4 describes the effects of a viscous layer on the visco-magnetic torques 
at the CMB and an Ekman number is estimated at the top of the core. The variations of the electrical 
conductivity in the lowermost mantle are studied in section 5. A final discussion ends the paper. 

2 MATHEMATICAL FORMULATION OF THE TORQUES 

At first order (©warning Citation 'poinlO' on page 4 undefined), the response of the rotating fluid 
core to Earth's nutations is a rigid body rotation. This approximation was checked experimentally 
(©warning Citation 'vany95' on page 4 undefined) and stays valid for large forcings (©warning Ci- 
tation 'noir03' on page 4 undefined). In the computation of the electromagnetic and viscous torques 
at the core mantle boundary (CMB), we may neglect the flow induced by the ellipticity of the CMB 
(©warning Citation 'sasa80' on page 4 undefined; ©warning Citation 'buff02' on page 4 undefined) 
and we describe the main flow in the outer core by an angular velocity ft. We consider the magneto- 
hydrodynamical equations in the frame of reference (e x , e y , e z ) rotating with the fluid outer core at 
the angular velocity vector fi, e z being defined by e z = Q/Q. The equations are made dimensionless 
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using as time scale, R the radius of the core as length scale and a typical magnitude of the radial 
component of the magnetic field B as magnetic field scale. The magnetic field and the flow velocity 
in the core (r < 1) are governed by the following dimensionless equations: 

— + (v-V)B = (B-V)v + £ m AB, (1) 
at 

Dv dft 

— + 2e z x v + — xr + fJx(f2xr) = -VP + E m A(V x B) x B + £Av , (2) 
where 

E 



qr 2 

is the Ekman number and v the kinematic viscosity. 

is the magnetic Ekman number which is Ekman number over the magnetic Prandtl number (y/rj) 
where rj is the magnetic diffusivity. 

is the Elsasser number, a = (por])^ 1 is the electrical conductivity of the core and p the density of the 
core fluid. 

In the above defined frame of coordinates , the motion of the mantle is a rigid body rotation Sum 
rotating at — e z defined by : 

Su> m (t) = 5lum [e x cos t — e y sin t] 

The angular velocity of the mantle is equatorial (no spin-up contribution) (©warning Citation 'buss68' 
on page 5 undefined; ©warning Citation 'noir03' on page 5 undefined). The dimensionless velocity 
in the mantle is described by: 

v M = Su M x r = -r5uj M [eesin(t + tp) + cos 9 cos(t + </>)] (3) 

where (e r , eg, e^) is the spherical coordinate system directly associated to (e x , e y , e z ). The magnetic 
field in the mantle (r > 1) is then described by the induction equation : 
<9B 

— + (vm • V)B = (B • V) v M + E% AB (4) 



where 



E M _ 



QR 2 

is the magnetic Ekman number of the mantle and tjm is the magnetic diffusivity of the mantle. 

When 5ujm = 0, the solution of (Q]) and dU is a diffusive poloidal magnetic field denoted Bo- As 
5u>m increases, magneto-viscous boundary layers develop around the core mantle boundary (r = 1). 
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The induced magnetic field in these boundary layers is denoted b. As 5u>m << 1, we have b « 
Bq = 0(1). The width of the magnetic skin layer at the top of the core (at the bottom of the mantle) 
is of order E m 1 / 2 ((E^ 1 ) 1 / 2 ) which is very small compared to 1. The viscous layer of size E 1 ! 2 is 
even smaller. Consequently, only radial derivatives of b and v have to be considered in the magneto 
hydrodynamic equations in the boundary layers. Moreover, we neglect the radial variations of Bo and 
vj\/ inside these thin boundary layers. A linearisation of equations (Q]),©>© with the above boundary 
layer assumptions leads to: 

dh d 2 <9v 

ot or z Or 

dv d 2 f dh dh \ 

Vr<l, — + 2e z xv-£-^v = -Vn + £ m A ((5— -B- — e r (6) 

ot or z \ or or J 

Vr>l, f-^£b = (7) 

where (3(9, (p) = B (r = 1) • e r is the radial component of the imposed magnetic field. 

Taking the curl of the motion equation © to eliminate pressure and using equation © to eliminate 
b, we obtain the following equation: 

,d ^d 2 ^d „ d 2 , ,„ . nl d „ d 2 ,dv „ ,„J(Vxv) 

Under our assumptions, V x v = — ^-ee + ^~ e <p an d it * s convenient to use the complex variables 
v + = vg + ivp, and v _ = vq — iv v to rewrite the last equation: 

d_ 

Or 

with the imposed velocity (in terms of complex variables) as boundary condition (r = 1) deduced 
from equation ©: 



d d 2 d d 2 d d 2 



E m AP 2 ^ , (8) 
or a 



vm± = — 2 — [(1 T cos @) ex P + _ (1 i cos 6*) exp — i(f + 99)] . 

The boundary condition imposes a time dependence of the form exp ±it which leads to a set of 
four differential equations deduced from ©. We use the exponent W to denote the sign of the time 
dependence and we define: 

v± = v±~^ exp it + v±^ exp —it 

In the following, we solve explicitly the problem for one component of the velocity v\_ ' and the other 
three components of the velocity may be easily deduced from it. The equation for v[ ^ is: 



D_ 

Or 



d 4 d 2 
EEm d^ + ( iE + iE ™ ~ 2iE ™ cos 6 ~ ^ A / ?2 )^2 + 2 cos 9 ~ 1 



v\ ) = . 



The associated polynomial function has four roots Z% corresponding to four exponential elementary 
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solutions which can be written (thanks to M. Greff): 

Zi = ±(2EE m )- 1 / 2 [-iE - E m {i - 2i cos 9 - A/3 2 ) 
± [-E 2 + 2EE m (l - 2cosfl - ik(3 2 ) + 



E 2 m (4 cos 9 - 1 - 2i A/3 2 - 4 cos 2 9 + 4i cos 0A/? 2 + A 2 /? 4 )] 1/2 



1/2 



As the velocity field must vanish far away from the boundary layer (lim rE -i/2_ ) ._ 00 v = 0), both 
roots with a positive real value are retained and the solution may be written: 

Vr < 1, v { ~ ] = X 1 exp(Z!(r - 1)) + X 2 exp(Z 2 (r - 1)) , 

where ATiand X 2 are constants to be determined. 
The equation (f5]) gives us the solution for b+ ' : 

Vr < 1, = - ^ X x exp(Z 1 (r - 1)) - X 2 ex P (Z 2 (r - 1)) . 

Using lim J .( B M)-i/2_ >QO 6 = 0, the solution for the magnetic field in the mantle may be deduced 
directly from ©: 

Vr>l, b { - ] =X 3 exp(Z 3 (r-l)) , 



where Z 3 = -(l + i)/y/2E%. 

We use the continuity of the velocity, the magnetic field and the electrical currents at the core 
mantle boundary (r = 1) to determine the constants X\, X2, X 3 . 

X!+X 2 = — ^-(l + cos0), 

" 0Zi (5Z 2 

-A 3 — ; 2 1 — ~ 2 2 j 

i E m Z^ i + E m Z2 

■ 2 1 — " 2 

i ~\~ E m Z^ % + E m Z2 



EmX 3 Z3 — E m 



With the solutions to this set of equations, the velocity and the magnetic field are fully determined 
within the boundary layers. With our scaling, the magnetic torque T m scales with pR 5 Q 2 E m A and the 
viscous torque T v with pR 5 Q 2 E. Here, we use the complex notation by introducing f = T x + iT y . 
The derivation of the viscous torque is given in Appendix A. 



T v = (e x + iBy) J J r x f v dS , (9) 

r f 2n dv-i. dv 

/ / [(l + cos#) — - + (l-cos9)— -}exp (iip) sin 9 d9 dip . (10) 
Jo Jo or or 



2 



The magnetic torque could be calculated by a surface integral (©warning Citation 'roch62' on 
page 7 undefined) similarly to the viscous torque. The magnetic torque could be deduced from the 
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perturbed magnetic field b at r = 1: 
f m = (e x + ie y )- [ rx (0b)dS , (11) 



P7T 



= - / /9[(1 + cos d)b+ + (1 -cos 0)6_] exp (i<p) sin fldfldyj . (12) 

It is of some use to introduce the coupling constant K deduced from the torque to compare with 
the observed data (©warning Citation 'math02' on page 8 undefined): 

« = , 

llOUJM 

where / is the dimensionless moment of inertia of the core. 

Both torques are integrated numerically using a (6, tp) grid where (3(6, tp) is prescribed. 

For a very weak magnetic field (the Lorentz forces tend to vanish) Toomre (©warning Citation 
'toom74' on page 8 undefined) predicted that the torque is pointing it /4 away from the direction of 
the imposed angular velocity Su. For very low Elsasser and Ekman numbers, the solution follows 
this asymptotic behavior (Irn(K) = —Re(K)). For a large and dipolar magnetic field, we compare 
successfully our results for a very low Ekman number (E = 10 -16 ) with the coupling constants found 
by Buffett (©warning Citation 'buff02' on page 8 undefined) in their inviscid study. In the limit of 
large Ekman numbers and low Elsasser numbers, we check that the computed torque tends toward the 
spin-over torque (©warning Citation 'gree68' on page 8 undefined). 



3 MAGNETIC FIELD AT THE CMB 

Even though the magnetic field at the CMB is dominated by the axial dipole component, all spherical 
harmonic components contribute to the electromagnetic torque. The spatial magnetic power spectrum 
(Mauersberger-Lowes spectrum) of the magnetic field at the Earth surface is deduced from obser- 
vations (©warning Citation 'lang82' on page 8 undefined; ©warning Citation 'voor02' on page 8 
undefined). At the CMB, this spectrum is fitted by the power law 1.085 10 10 (0.959)'nT 2 (©warning 
Citation 'stac92' on page 8 undefined) if one excepts the dipole component (I = 1). This is a rela- 
tively flat magnetic spectrum which has to become stiffer at very large I to satisfy energetic arguments 
(©warning Citation 'robe03' on page 8 undefined; ©warning Citation 'chri04' on page 8 undefined). 
Taking the same electrical conductivity for the core and the mantle (E m = E% = 1.8 10~ 9 ) and 
a very low Ekman number (E = 10 -16 ), we compute the electromagnetic coupling constant using 
the first l max degrees of the magnetic field (higher degree coeficients are set to zero). Figure Q] illus- 
trates possible contributions of the small scales (I > 13) of the magnetic field to the amplitude of the 
electromagnetic torque. Different random sets of spherical harmonic coefficients matching the spectra 
dependence, give comparable contribution to the torque (variations lower than 10%). Consequently, 
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max 

Figure 1. Imaginary part of the electromagnetic coupling constant versus the truncature level l m ax of the mag- 
netic field at the core boundary. The observed geomagnetic spectra is extrapolated randomly with a power law 
(0.959)' for I > 13 and the coupling constant is computed with the visco-magnetic model (VMM) or with the 
weak field model (WFM). Each line represents a set of spherical harmonic coefficients sastifying the spectra 
dependence. The bold lines show the mean values of both model. 

the mean value of the coupling constant is representative of what could happen at the CMB and in the 
following, we keep only the mean value to present the results. 

The contribution of the large degrees of the spherical harmonics of the magnetic field in the visco- 
magnetic model is smaller that the one associated with the weak field model (©warning Citation 
'buff92' on page 9 undefined). The back reaction of the Lorentz forces on the flow is to reduce the 
electromagnetic torque at the boundary. This effect is emphasised for the small scales of the magnetic 
field. In some cases (for example, all coefficients positive), the contribution of the large degree is 
negligeable and the coupling constant curve becomes flat (highest curve in Figure [T]). 

Figure [2] shows that the mean electromagnetic torque associated with the observed magnetic 
field at the CMB ((0.959)') is too low to fit the imaginary part of the observed coupling constant 
(—1.85 10~ 5 ). Following the ideas of Buffett (©warning Citation 'buff92' on page 9 undefined; 
©warning Citation 'buff02' on page 9 undefined), we explore the effect of an increase of the magnetic 
field at small scales . We studied two different power laws for the spectra (I > 13): (1.044)' giving an 
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Figure 2. Imaginary part of a mean electromagnetic coupling constant versus the truncature level l max of the 
magnetic field at the core boundary. The coupling constant is computed with the visco-magnetic model with 
three different extrapolated spectra for / > 13. 

magnetic energy 10 times greater than the standard one at degree 40, and (1.098)' corresponding to an 
energy 10 times greater than the standard one at degree 30. An increase by a factor 10 in energy at de- 
gree 40, increases the overall coupling constant less than a factor 2 which remains too small to match 
the observed value. In Figure |2l we see that only the (1.098)' spectrum could explain the observed 
data. From a geophysical point of view, this spectrum is unlikely as it dissipates a large amount of en- 
ergy. Using the result of Roberts et al. (©warning Citation 'robe03' on page 9 undefined) (eq 2.7 page 
104), we found a ohmic dissipation of 0.037W for I < 40. This is large compared to the dissipation 
associated to dipolar component alone which is 0.08GP^. According to the scaling deduced from nu- 
merical dynamos (©warning Citation 'buff02b' on page 9 undefined; ©warning Citation 'robe03' on 
page 9 undefined; ©warning Citation 'chri04' on page 9 undefined) which takes into account the dis- 
sipation of the toroidal part of the magnetic field and the contribution of small scales of the magnetic 
field, this (1.098)' spectrum is too dissipative. For example, with a dissipation of 1.32GW associated 
to the large scale magnetic spectrum, Roberts at al. (©warning Citation 'robe03' on page 9 undefined) 
estimate a total ohmic power loss between 1 and 2TW. 
Another dissipative process is thus needed to explain the nutation data. 
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4 VISCOUS EFFECTS AT THE CMB 

In this section, we assume the electrical conductivity in the core and in the mantle to be the same 
(©warning Citation 'buff02' on page 11 undefined) and focus on viscous effects. Figure [3] shows 
components of the velocity and perturbed magnetic fields in the boundary layers for E m = E%[ = 
1.8 1CT 9 (see tableQ]) and different Ekman numbers. For very low Ekman numbers, the viscous layer 
is very narrow (E 1 / 2 ) and the magnetic field (symmetry and amplitude) is nearly unchanged by the 
presence of the Ekman layer. For Ekman numbers comparable to the magnetic Ekman number, the 

1/2 

width of the viscous layer becomes as large as the magnetic skin depth (Em ) and magnetic field is 
induced deeper into the core. Consequently, the perturbed magnetic field looses its symmetry and its 
value at the CMB decreases. 

This physical behavior is summarised on figure [4] The magnetic torque does not vary for very 
low Ekman number (E < 1CP 11 ) and both components decrease as the Ekman number approaches 
the magnetic Ekman number. As expected, the viscous torque increases with the Ekman number. The 
imaginary part of the magnetic and viscous coupling constants become comparable for E 2. 1CP 12 
while their real parts match for a larger Ekman number (E 4. 1CT 10 ). This difference results directly 
from the geometry of the spin over viscous torque which exhibits a very low imaginary part (0.259) 
compared to the real one (2.62) (©warning Citation 'gree68' on page 1 1 undefined). As a conclusion, 
the resulting torque at the CMB is largely modified by a viscous shear layer for E > 1CT 12 even if 
the depth of the viscous layer is much smaller than the magnetic one (Figure [3^) and its effect barely 
changes the induced magnetic field at the boundary (Figure 03). 

In order to fit the imaginary part of the observed coupling constant (—1.85 10~ 5 ), an Ekman 
number of 3. 10~ n is necessary as shown on figure 01 With such a value, the viscous dissipative 
process represents 85% of the whole dissipation at the CMB but the real part of the coupling constant 
is still dominated (75%) by the magnetic torque. 

5 ELECTRICAL CONDUCTIVITY AT THE BOTTOM OF THE MANTLE 

In the visco-magnetic model for CMB parameters, the influence of the electrical conductivity of the 
lowermost mantle is secondary. Figure [5] shows the evolution of the coupling constant as the electrical 
resistivity is increased up to 10000 times the electrical resistivity of the core for an Ekman number 
of 4. 10 . The imaginary part of the coupling constant varies less than 20% but the real part of the 
coupling constant is divided by 3. As expected, for large magnetic Ekman number in the mantle, the 
torque is mainly dominated by the viscous part of the torque. 

A trade off between viscous and magnetic torque could be found in order to fit the observational 



12 B. Deleplace and P. Cardin 



0.00 



-0.05 



^ + -0.10 
5 -0.15 



-0.20 



-0.25 



E=10" 

— E=10" 11 

E=10~ 9 

E=10" 8 



-0.30 
-0.0004 



-0.0003 



1 



V 

/ \ 



v 

A 



-0.0002 -0.0001 
radius (r-1) 



0.0000 




Figure 3. Velocity Im(v ( + ) and magnetic | b\_ ' | perturbed fields at the core-mantle boundary for E, 
1.8 10~ 9 and different Ekman numbers. 
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Figure 4. Variations of the coupling constant K as a function of the Ekman number. K v (K m ) is the viscous 
(magnetic) component of the coupling constant 



data of nutations. On figure [6l for each value of the electrical conductivity of the mantle (E%[), we 
plot the Ekman number (E) corresponding to a total torque in agreement with the observational con- 
straint Im(K CMB ) = —1.85 10 -5 . For the nearly flat standard spectra (0.959)' (corresponding to the 
magnetic field at the CMB), Ekman numbers between 2 and 5. 1CT 11 are retrieved from the inversion 
whatever the conductivity at the bottom layer of the mantle. For the largest increasing spectra (1.098)', 
Ekman numbers vary more significantly with the conductivity of the lowermost mantle and very low 
Ekman numbers are retrieved when the electromagnetic torque becomes significant (comparable elec- 
trical conductivity on both sides of the CMB). 

The electrical conductivity at the bottom of the mantle is difficult to determine. Theoretical anal- 
ysis and experimental measurements indicate that silicate rocks have a lower electrical conductivity 
than the liquid metal of the core (©warning Citation 'poir92' on page 13 undefined; ©warning Cita- 
tion 'shan93' on page 13 undefined). Discoveries of new phases of perovskite, such as post perovskite 
(©warning Citation 'iita04' on page 13 undefined), or metal alloys of silicates may change this state- 
ment. The resistivity may even present lateral variations as shown by seismic lateral variations in the 
lowermost mantle (©warning Citation 'lay98' on page 13 undefined). Except in the case of both a 
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Figure 5. Variations of the coupling constant with the Ekman magnetic number of the mantle. 

large electrical conductivity at the base of the mantle and a large magnitude of a hidden magnetic field 
at the CMB, the Ekman number needed to explain the observational data is around 10 -11 . 

6 DISCUSSION 

The real part of the coupling constant is not dissipative and directly influences the period of the nu- 
tations (©warning Citation 'deha97' on page 14 undefined; ©warning Citation 'hindOO' on page 14 
undefined; ©warning Citation 'math02' on page 14 undefined). But the discrepancy between the the- 
oretical and observed periods of the free core nutation (FCN) is too large to be explained only by the 
real part of the coupling constant associated with the dissipative torque at the CMB. Hence, a dynamic 
ellipticity of a few hundred meters at the CMB has been introduced to account for this discrepancy. 
Then, the real part of the coupling constant cannot be used anymore as an observational constraint 
to determine the nature of the dissipative torque at the CMB. It is true, though, that a visco-magnetic 
dynamic model of the CMB reduces the real part of the coupling constant compared to a weak field 
model and consequently tends to increase by 10% the estimate of the dynamic ellipticity at the CMB. 

The visco-magnetic model of the magnetic skin layer shows that the small scales of the magnetic 
field at the CMB, that cannot be directly inferred from magnetic observations, do not contribute sig- 
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Figure 6. Curves in the plan (E, E™ ) on which Im(K CMB ) = -1.8510~ 5 , for E m = 1.8 10~ 9 and the three 
different spectra for (I > 13). The observed magnetic field at the CMB corresponds to the nearly flat spectra 
(0.959)'. 

nificantly to the electromagnetic torque between the core and the mantle. With a fiat or decreasing 
power magnetic spectrum at the CMB, the electromagnetic torque is too weak to explain the coupling 
constant Im(K CMB ) = —1.85 10 -5 , even if the electrical conductivity of the lower most mantle is 
comparable to the core one. 

The visco-magnetic model of the boundary layers at the CMB using the observed geomagnetic 
field imposes the presence of viscous dissipative effects. Apparent Ekman numbers between 2 and 
4 10~ n are needed to fit the observational constraint (corresponding to a turbulent viscosity of 
3.5 10~ 2 m 2 s~ 1 ). This observational constraint may change as the quality of the data and their treat- 
ment improve (©warning Citation 'florOO' on page 15 undefined). Only a reduction of the value of 
Im(K CMB ) by a factor 3 would make viscous effects unnecessary to explain the observations. How- 
ever, the recent study of Palmer & Smylie (©warning Citation 'palm05' on page 15 undefined), if 
correct, gives an Ekman number of 7 10~ n using their own analysis of VLBI data of the free core 
nutation and a pure viscous model of coupling at the CMB. The agreement is also very good with the 
results of Mathews & Guo (©warning Citation 'math05' on page 15 undefined) which states that an 
Ekman number larger than 5 1CT 11 is needed at the CMB. Such values of the Ekman number at the 
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CMB are also compatible with the dissipation needed to account for the relaxation of torsional oscil- 
lations in the core (©warning Citation 'zatm97' on page 15 undefined; ©warning Citation 'jaul03' 
on page 15 undefined). 

At this stage, we like to see this observational constraint as a measurement of an apparent viscosity 
at the top of the core. The visco-magnetic model suggests an effective viscosity five thousand times 
larger than the expected molecular viscosity of iron at the core conditions. We would like to stress than 
an apparent viscosity is space and time dependent whereas, it is defined here on the whole surface of 
the CMB at the diurnal frequency. Consequently, this value of viscosity (3.5 10~ 2 m 2 s~ 1 ) may not 
be generalized to the bulk of the core, at small scales and at different time scales. This is large but 
comparable to effective viscosities used in fluid dynamics of the ocean or the atmosphere. In these 
fields, a large apparent viscosity is the net result of the turbulent transport due to the small scales of 
the flow on the large scale flow. Such an explanation may be valid in the Earth's core even though 
we do not have any evidence for the action of small scales at the CMB. As discussed by Davies & 
Whaler (©warning Citation 'davi97' on page 16 undefined), the effective transport of momentum 
may be generated by convective motions associated to the dynamo process or by surfacic flows such 
as topographic winds, unstable boundary layers motions, or chemical/compositional fluxes. 

This work has been financed by the program DyETI and PNP of CNRS/INSU. The authors would 
like to thank Thierry Alboussiere, Dominique Jault, Marianne Greff and Veronique Dehant and the 
rewievers for very helpful comments and suggestions. 

APPENDIX A: TORQUE FORMULATION 

The viscous torque is computed from the viscous forces at the core mantle boundary: 



where i v is the viscous force per unit area. In our geometry and within the boundary layer approach, 
the viscous force on a sphere may be written : f v = f Vrf) eg + /„ e v = ^p-eg + ^re^. At r = 1, we 
have: 

r x f„ = {-f Vrip cos 9 cos <p - f Vr9 sin (p)e x + (-f Vrip cos 9 sin (p + f Vr0 cos ip)e y + f Vrip sin 9 e z 
Introducing f = T x + iT y , the dimensionless complex viscous torque at the CMB is: 





which in terms of v + and V- may be expressed: 




+ (1 — cos 9) — — ] exp (iip) sin 9d9d(p 
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which is the expression shown in equation (fTOb . Equation (fTTb being similar to equation Q, the deriva- 
tion of the magnetic torque is similar to that of the viscous torque. 



